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Abstract 

We investigate behavior of a noncritical model in external electric field 
and explore its phase structure in the quenched approximation Nj <^ N^. 
We compute the conductivity of QCD plasma in this model and compare it 
with the predictions of Sakai-Sugimoto model, D3-D7 system and the lattice 
simulation. We find that, while the behavior of conductivity in noncritical 
model as a function of temperature and baryon density is similar to those of 
D3-D7 system, the phase diagram of noncritical model resembles the phase 
diagram of Sakai-Sugimoto model. 



1 introduction 



The AdS/CFT correspondence is a duality between the strongly coupled conformal 
field theories and string theory in a higher-dimensional Anti-de-Sitter space-time [T]. 
The generalization of AdS / CFT correspondence to more realistic gauge theories like 
QCD, provides new insights to understanding the dynamical non-perturbative effects 
in QCD, such as confinement, chiral symmetry breaking, color superconductivity and 
so on. 

Dual gauge theories arising from brane constructions in ten-dimensional criti- 
cal string theory are supersymmetric. In order to break supersymmetry, one may 
compactify the supersymmetric guage theory on a circle of radius R and impose 
anti-periodic boundary conditions for fermions around the circle. The resulting ef- 
fective theory at low energy compared to the Kaluza-Klein mass scale, Mkk ~ ^, 
is pure QCD without fundamental matter. Fundamental matter may be incorpo- 
rated by adding fiavor branes [2]. In principle, it is possible to extract the physics 
of strongly correlated QCD by using these holographic models. The main obstacle 
of this approach, however, is that resulting holographic QCD contains undesired 
Kaluza-Klein modes with the mass as the same order of hadrons and glueballs. One 
way to overcome this problem is to consider brane backgrounds in noncritical string 
theory [13]. Since in this case holographic backgrounds live in lower dimensions, the 
problem of extra KK modes is more tractable . 

A noncritical holographic model of QCD has been introduced in [3], where a 
stack of Nc D4 branes in six-dimensional noncritical string theory play the role of 
color branes. Fundamental matters may be added by inserting Nf D4 branes and 
D4 breanes in this background [1] . Further investigations [HIS] have shown that this 
model captures correctly many properties of low energy QCD, like mass spectrum 
of mesons, area low behavior for the Wilson loop, glueball mass spectra, etc. In par- 
ticular, a holographic calculation of Wilson line has been performed recently in [H] 
for near extremal D3 branes, D4 branes, non-critical near extremal AdSe model and 
Klebanov-Strassler model. It has been shown that the noncritical background ad- 
mits a reasonable fit to lattice results (Actually the Klebanov-Strassler background 
exhibits the best fit to lattice results, but noncritical model has a asymptotic AdSe 
metric, see for details [E]). The thermal phases of the model have been investigated 
in [6] , where it was shown that there is a first order confined / deconfined phase tran- 
sition at Tc = In similar to Sakai-Sugimoto model, chiral symmetry restores 
above a critical temperature T^sb = ^^^i^, where L is the separation of flavors branes 
at infinity. 

It is also interesting to investigate the behavior of QCD matter in presence of 
other fundamental interactions. Indeed the physics of some astrophysical phenomena 
like quark stars and notron stars is related to the behavior of thermal QCD in 
external electromagnetic field. For example conductivity of QCD matter encodes 
the timescale for expulsion of magnetic fiux lines from the core of quark star (see for 
example [I5]). Since QCD is strongly coupled at these circumstance, it is natural 
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to use AdS/CFT and extracting the properties of QCD at these high densities and 
temperatures. Using gauge/garvity duahty, response of Sakai-Sugimoto model to 
external electric and magnetic field has been studied in [7]-[l0]. behavior of massive 
M = 2 hypermultiplets in an TV = 4 SYM plasma (D3-D7 system) in presence of 
background electric and magnetic field was explored in [HKHKin] and the effect of 
magnetic field on dynamics of noncritical holographic model has been considered 
in [17]. 

Motivated by these discussions, in this paper, we consider a noncritical holo- 
graphic model of QCD, and examine its behavior in the presence of external electric 
field. We compare the results of this model with the predictions of Sakai-Sugimoto 
model [7], D3-D7 system [16], and a lattice simulation [19]. We find that the be- 
havior of conductivity in noncritical model as a function of temperature and baryon 
density is similar to the results of D3-D7 system. In weak field regime, where lattice 
simulation is accessible, noncritical holographic model predicts a finite conductiv- 
ity for deconfined QCD matter which is linear in temperature, in good agreement 
with lattice QCD result. Also, the phase diagram of noncritical model resembles 
the phase diagram of Sakai-Sugimoto model. In particular we find that electric field 
reduces chiral-symmetry restoration temperature. 

This paper is organized as follows: In section 2 we briefly review the noncritical 
model. In section 3 we analyze the dynamics of flavor branes in presence of electric 
field and extract the phase diagram of the model. We check our result by making 
use of Kubo formula in section 4 and section 5 is devoted to brief summary and 
conclusions. 

2 Review of the model 

In this section we briefly review the noncritical model of [3]. The model is based on 
D4/D4-D4 brane system, where N^. D4-branes compactified on 5*^ with radius R and 
Nf D4-D4 flavor branes are transverse to the S^. By imposing periodic boundary 
condition on bosonic field and antiperiodic boundary condition on fermionic field 
along the circle 5*^, supersymmetry is broken and one obtains QCD spectrum at the 
energies below the Kaluza-Klein scale. The brane configuration of the system is 

t Xi X2 X3 X4 X5 

D4 X X X X X 
D4-D4 X X X X X 

where X4 is the coordinate of S^. At zero temperature the 6-dimensional back- 
ground metric is given by 
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dsl=(-^)\-dt' + dx'dx' + f{u)dxl)+(^y^ (2.1) 

\R-AdS / V ^ / JK^J 

/H = l-(^)^ RAds=^is 

There is also a constant dilaton and a 6-form field strength 



f u Y 

F(6) = —Qc I ^ I dxo A (ixi A dx2 A (ixa A (ix4 A du 

The space in (0:4, u) plane looks like a cigar where tip located at n = ua, and to 
avoid a conical singularity the periodicity of x^^ should satisfies 

X4 ~ X4 + 2nR = X4 + ^^^^^ (2.2) 

Kaluza-Klein energy scale is Mkk = -5 = gp""^ , below this scale the effective 
theory is QCD4. At nonzero temperature there are two solutions with the same 
boundary condition [6]. The one which dominates at low temperature is given by 
(12. ip where the periodicity of euclidian time is arbitrary, tE ~ tE + f3, and periodicity 
of X4 is given by (12.21) . The flavor branes form a U embedding configuration, so this 
background corresponds to the confined phase with broken chiral symmetry. 

By increasing temperature a confinement / deconfinement phase transition occurs 
at Tc = and the background for Tc <T is represented by 

d-^' = {-^fUWtl + dxdx^ + dxl) + (2.3) 

RAdS U J [Uj 

f[u) = l-[ j, lE^lE-l 1 

U ' OUT 

In this case the periodicity of X4 circle is arbitrary and temperature is given by 
T = . J . This background allows two embeddings for D4-D4 flavor branes: U 
embedding which is preferred configuration for < T < T^sb and parallel embed- 
ding which dominates for T > T^5b, where T^sb = Thus chiral symmetry 
restores at temperatures above T^sb- 



3 Adding U{1) gauge field 

In order to accommodate an external electric field, following [16], we turn on a t/(l) 
gauge field on flavor branes 
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A^{tE,u) = -iEtE + h{u) 



(3.1) 



where time dependence part of gauge field describes a static electric field on the 
boundary and u dependence part encodes the response current. In this section we 
analyze various phases of noncritical holographic QCD in external electric field by 
making use the method of [16] and [7]. 



3.1 Deconfined Phase 

Let us start with investigating the effect of an external electric field on high tem- 
perature (T > Tc) phase of QCD. By using ^ = {tE,^,u) for parametrization of 
D4-D4 flavor branes, the induced metric on flavor branes is given by 







ds' = {^)\f{u)dtl + dxdx.)+ (i^fx',' + {^Y-^)du' (3.2) 

RacIS V RAdS U f{u)J 

DBI action in presence of background gauge field (13. ip takes the following form 



(3.4) 

where M = 2N fT^e'"^ , e = 2TTa'E, and a'^ = 27ra'A^. By making use of the 
first integrals of motion one can find the asymptotic forms of gauge field and x'^ at 
large u 

< ~ ^, 2:4 ~ ^ (3.5) 

(3.6) 

where j and c are identified with the current and chiral condensate in dual gauge 
theory. By doing holographic renormalization one can show that the physical current 
on the boundary is 

^' - (3.7) 

^AdS 

By writing the action in terms of current, j, we have 



5 = ^/rf^«y(/(«)xf + ^)(/(«)-^)(/(«)-^)"i 



^Throughout this paper we work with action density 
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As we have mentioned in the previous section, deconfined background possesses 
two embeddings, U embedding and parallel embedding. Our aim is to determine 
current j for a given external electric field in these embeddings and then extracting 
conductivity of these configurations. 

First consider the U embedding. In the case of vanishing response current, j = 0, 
the solution takes the following form 



2:4 (m) 



^AdS 



u- 



2 p4 

-ft / 



- 1 



(3.9) 



where Mq is the turning point, ^4(^0) 
action we arrive at 



0. By inserting the above solution into the 



u 



Af 

^AdS Juo 



du U 



^ ^AdS 



u^f{u) 



1 - 



g -"-AdS ' 
^0 ' 



(3.10) 



As long as Cg < -^r—u^f^uo), the action is real and embedding is acceptable. 
Since by turning on the current, action increased, the fevered configuration is a U 
embedding with vanishing current, j = [7]. In dual QCD this means that the 
deconfined chiral-symmetry breaking phase is an insulator. 

A natural question is that what happens for e > cq?. Figure 1 depicts L as a 
function of c for different values of electric field. As it is evident from this figure, L 
is a decreasing function of e, so there is a maximal value of e at fixed values of L and 
T, such that above which there are no U embedding solution. Thus we expect that 
the favorite solution in this regime becomes parallel embedding and there should be 
a phase transition from chiral-broken phase to chiral symmetric phase by increasing 
electric field [7]. We come back to this issue after discussing the behavior of parallel 
embeding in electric field. 

In the parallel embedding the DBI action becomes 



Su 



^AdS Jut 



du U 



fin) 



p2 d4 



(3.11) 



From this expression, it is clear that action becomes complex somewhere unless 
a nonzero current being turned on. The magnitude of current is given by 



3 — ^AdS ^ 



where Uc is the root of numerator 

f{Uc 



^ ^AdS 



(3.12) 



(3.13) 
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Figure 1 : L as a function of c in deconfined phase for different values of electric field 

Therefore the chiral-symmetry restoration phase behaves like a conductor by 
following conductivity 

_r _ (27raOW j _ (2W)W 

^ non— critical j-, j-,3 j-, ^cl,*^) J v'^'-'-^/ 

There is no algebraic solution for f l3.13p . however, we can study the weak filed 
and strong field behavior of conductivity 



'i|^(27ra')^Ar T E<^^T^ 

^non— critical \ (3.15) 

It is interesting to compare this results with the predictions of Sakai-Sugimoto 
model, (Js-s, and D3-D7 system, aos-m- The weak and high field behavior of 
conductivity in these models are given by [TtlTB] 

T2 E « ^AsT^ ( « ^v^T^ 



where A5 and A4 stands for five and four dimensional t' Hooft coupling. As 
it is clear from above expressions, the behaviour of conductivity as a function of 
temperature and electric field in noncritical model is similar to the D3-D7 system. 

Now, let us consider to the phase structure of deconfined background. Since 
external electric field bring the system out of equilibrium, as argued in [7], to de- 
termine which phase is favorite, one can use a Maxwell Like construction for order 
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parameter c = ^ |e,T- According to Maxwell construction, for a fixed value of L 
and T, the transition occurs at the value of e such that two areas A and B (see 
figure 2(a)) become equal. By changing electric field we find the phase diagram 
as figure 2(b), which is, remarkably, similar to those of Sakai-Sugimoto model [7]. 
According to this phase diagram, we observe that the critical temperature decreases 
with increasing external electric field, as we expected from the polarization effect of 
electric field. Also at zero electric field we get the result of [6] for chiral symmetry 
restoration temperature, T^b = 



L 




t 0.00 I ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ^ ' 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 0.2 0.4 0.6 0.8 1.0 



Figure 2: (a) Maxwell construction (b) phase diagram fort L = 1 in presence of 
external electric field in deconfined background 



3.1.1 finite density 

In this subsection, we generalize our analysis in the parallel embedding for finite 
baryon density by turning on a nontrivial zero-component of the gauge field ao{u). 
The resulting DBI action has the form 



Solving the equation of motion one finds the asymptotic behavior of gauge fields 

as 

j d 
Qx — constant qq constant ^ (3.16) 

where j and d are related to the physical current and charge density via 

{2na')M . (W)Ar 
J = —3 J D = —3 d (3.17) 

-^AdS -^AdS 

Writing the action in terms of current and charge density we arrive at 
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s 



^AdS 



duw"' 



\ 



p2 p4 



p-d^f{n) 



Reality of action implies a nonzero current 



from which conductivity at finite charge density can be read as 



(3.18) 



(3.19) 



J_ _ {2Tia'fU 

E RAdS 

At high density approximation it becomes 



rf2 



'U 



c ' 4 



(3.20) 



_ 25a^W d _ 5 D 

'^noncntical 4J^5 127r (^O.ZiJ 

Let us compare this result with that one obtained in Sakai-Sugimoto model and 
D3-D7 system 



_ 27 D _ 2 D 

(^Sakai-Sugimoto — — 777—-^ O'DS-D? — J=-^ (3.22) 

87r^A5 T-^ 7rvA4 

Again we observe that charge density dependence of conductivity in noncritical 
model and D3-D7 system is the same. 



3.2 confined phase 

In this section we will study the response of confined phase to external electric field. 
This background dominates for T < T^. and in the case of absence external field, 
the only allowed embedding for flavor branes is U embedding . To proceed we start 
with the DBI action for flavor branes 



S=^l duu'Jif{u)x',' + %)(1 - ^) + (3.23) 

^AdS J ^ U U 

First consider the U embedding with j = 0. The equations of motions gives 



^AdS 

u^fiu) 



r^io 



f{u){l 



</(^o)(l- 



p2 d4 



- 1 



AdS 



(3.24) 



8 



with the following asymptotic behavior near the boundary 



(3.25) 



where 



p2D4 



(3.26) 



In order to have a physical solution, it is required mq > e^RAds- Since the range 
of coordinate u is (ukk, oo), if e < -nf^, uq can interpolates between ukk and oo. 
Figure 3(a) shows the behavior of L in terms of c, in this case. On the other hand, 
for e > -^f^, the range of uq would be (e^RA^SyOo) and the plot of L(c) is like 

2 

deconfind background (figure 3(b)). Then for e > , there is a maximal L for a 
given e. Since by increasing e the maximal of L decreases, there is a maximal e at 
any fixed value of L, such that for e > Cmax, U embedding is not a valid solution. 
As discussed in [7] the solution for e > Cmax is a V shape embedding, where flavor 
branes follow parallel radial geodesies. Thus this solution satisfies 0:4 (n) = except 
at the tip. DBI action for D4-D4 branes in this background is as follows 





2 2 

Figure 3: L as a function of c in confined phase (a) e < (b) e > -^f^ 

^ ' ^AdS ^AdS 



^AdS 



U 



du ■ 



^ ^AdS 



I - 11 



From which one can read the conductivity of V embedding 

{2na')^Af E'2 



a 



(3.27) 



(3.28) 



Therefore the confined phase exhibits two phases in external electric field. U 
embedding which describes an insulating phase and V embedding with a finite elec- 
tric conductivity. In order to explore the phase diagram, we apply the same method 
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used for deconfined phase. The resulting phase diagrams for fixed value of ukk and 
L is shown in Figure 4, which are similar to those of Sakai-Sugimoto model [7]. From 
this diagram we observe that U embedding is dominate background for e < pf ^ , 
disregard how much L is, as we expected. 



20 r 4r 




0.00 0.02 004 0.06 0.08 0.10 0.0 0.5 1.0 1.5 



Figure 4: phase diagram in confined phase, (left) ukk = -5, (right) L=l 



3.3 electric susceptibility 

The response of an insulating phase to an external electric field can be measured by 
electric susceptibility, which defined by 

Xe = (3.29) 

Since there is an insulating phase in both deconfined and confined background, it 
is interesting to use holographic description to calculate the susceptibility in the dual 
QCD. Indeed holographic free energy is divergent and one has to add a counterterm 
to find finite results. We apply regularization of [7] 

Xe = + ^ le^o (3.30) 

Using this expression, we find electric susceptibility in deconfined and confined 
phase according to figure 5. 



4 Kubo formula 

The response of a thermodynamic system to an applied external field is described 
by transport coefficients of the system. For small deviation from equilibrium, Kubo 
formula relates transport coefficients to the equilibrium retarded green's functions 
of the system. In particular real-time correlator of two electromagnetism current 
determines electric conductivity of the medium via 
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0.00 0.01 0.02 0.03 0.04 0.05 0.00 0.01 0.02 0.03 0.04 0.05 0.06 



Figure 5: susceptibility in deconfined phase (left) and confined phase (right) 



a= \ hm i/mG^^(a;,k) (4.1) 

where Gaa is retarded green's function of two transverse electromagnetic current. 
In this section we calculate retarded green's function of two EM currents in the dual 
QCD by using of the Lorenzian AdS/CFT prescription [TH|. To do so, we turn 
off external electric field, e = 0, thus the system is in equilibrium and turn on a 
nontrivial time- component gauge field Aq as a background field, corresponding to a 
finite charge density in the dual gauge theory. Therefore the background is given by 



ds^\r>. =(^)\-fiu)dt' + dx^dx^)+(^y^ (4.2) 

\RAdsJ \ U J f{u) 

Ftu =-A[{t,u) 

Where ds^lm is the induced metric on D4-D4 flavor branes. According to holo- 
graphic dictionary we focus on the linearized fluctuations of U(l) gauge field on the 
gravitational background. By expanding the DBI action up to second order in the 
field strength around the background (14.21) . we arrive at 



(W)W [ du . ju' + C') i^^' + C'Y , 

^Rl,s J V^^^TC-^^ fin) ^ u > 



where, C = ^^^^^^ D, it is convenient to change variable to 



u 
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(14. 3 p take the following form 



2 J Jl + C'^y^ ^ f{y) t1 ' 



+/(y)(i+c^'V)E4] 



where vt = -^^"^ = 7% and 



C'^^ (4.5) 

Since transverse part of electromagnetic current couples to transverse part of 
bulk electric field, we consider the e.o.m for transverse part of electric field. It is 
convenient to work in fourier space 

Am{x^ y) = j -0^/'-'' AM{k,, y) k, = (-a;, g, 0, 0) (4.6) 

and equation of motion for transverse components, Ea = ujA^ {a = 2,3), be- 
comes 



2 Q^f(u) 

Ul — 



K + , J^—^ iy-'mVTTC^yE'^ + ^^^E^ = (4.7) 

/(?/)v/i + c"V f(y) 

From this equation the near horizon behavior of E^ can be read as 

E. = ivT - y)^™/' (4.8) 

where w = = and ± represent ingoing and outgoing wave into horizon, 
respectively. By imposing ingoing boundary condition at the horizon, the behavior 
of solution near y = becomes 

E^ = A + By^ (4.9) 
The relevant part of boundary action for calculating retarded corellator is 



(4.10) 



(2W)W„ f dwdq 2AB 

R-AdS / 77r~v? — 9~ (4- 11) 



2 ""^7 (27r)2 u;2 
By applying AdS/CFT reception for calculating real-time correlator, we have 
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= j£k. = ^'lik. = 2(2™')WiS.../™|^/B| (4.12) 



In appendix we have calculated A/B at low frequency, w — > 0, the result is 
substituting this result into fl4.ip we arrive at 



l + _,^/l + C"y6_ (4.13) 




cr = ^— ^ UtJI + ^ (4.14) 



which is exactly f l3.20p with e = 0. 

5 Conclusion 

We have studied a noncritical holographic model in external electric field and com- 
pare its results with the ones of Sakai-Sugimoto model and D3-D7 system. We found 
that, remarkably, the behavior of conductivity as a function of temperature, electric 
field and baryon density is similar to the result of D3-D7 system. In particular in 
weak-field regime the conductivity grows linear with temperature, a oc T. 

Computation of transport coefficients by using lattice QCD requires an ana- 
lytic continuation to real-time space, which leads to systematic errors in results. 
However, lattice results with small errors provide insights about strongly coupled 
regime of QCD. In [19] putting systematic errors under control, a lattice simulation 
of conductivity in QCD has been done, the reported result is 



.rr.. I 7.5 ±0.8, T = 1.5T, 

' -CemI 7.7 ±0.6, T = 2% (5.1) 



T 



7.0 ±0.4, r = 3r, 



c 



where electromagnetic vertex factor is given by Cem = 47ra which a is 

fine structure constant and e/ is electric charge of a quark with fiavor f {Cem ~ ^ 
for two fiavors [19]). If we can trust this result, it shows that, with good accuracy, 
conductivity is linear in temperature, in accordance with the predictions of noncrit- 
ical models and D3-D7 system (soft wall model also predicts a linear dependence 
a ocT [2I]). 

By studding the phase structure of the model, we observed that the electric field 
reduces chiral-restoration temperature. Also the general structure of phase diagram 
closely resemble phase diagram of Sakai-Sugimoto model [7]. In addition we have 
checked our results by using Kubo formula in section 3. 
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Appendix 

In this appendix we explore the low-frequency behavior of retarded green's function, 
by solving fl4.7p perturbatively in uj. Let us first rewrite f l4.7p in terms of more 
convenient variable 

E^ = {yT-yr^/^{y) (5.2) 
substituting this into (14. 7p we find the following equation for 

0"(y) + ^0'(y)+fH0(y)=O (5.3) 

where we have rescaled variables as y — y yxy and C — 5- C = y^C, for 
simplicity and 

^ =J^+ (1 + 4?/5 - 2Cy6 + rcyii) 



l-y i-l + y)y{l + y + y^ + y^ + y^)il + Cy^) 
w{-2i + w) 25w'^y^{l + Cy) iw{l + Ay^ - 2Cy^ + 7Cy^^] 



4(-l + ?/)2 4(_i +^5)2(1 +c'?/6) 2{-l + y)^y{l + y + + + y^){l + Cy^ 

In order to solve this equation perturbatively in w, we consider following expan- 
sion for 

(j) = (j)o + uj(j)i (5.4) 

plugging this into (15.31) and imposing regularity condition at the horizon, we find 
that 00 is a constant and (pi satisfies following equation 

cP'liy) + Piy)<P[iy) = Riy) (5.5) 

where 

Piy) .1±V^%!£±Z^ (5.6) 
{-y + y'){l + y^C) 

Riy) -.[--^ + l + %^-2.-C + 7,nC 

2(-l + yy 2(-l + y)2y(l + y + y2 + y3 + ^4)(1 + ,^6^) 



14 



The general solution of above equation is 

My) = A0(y) + / dx^ - / dx 'J;] (5.8) 

Jo cj)'[x) Jo (f)'{x) 

where A is a constant and particular solution, (j){y), is given by 

0(y) = r dxe- ^0 '^'^^'^ (5.9) 



For our purpose, it is sufficient to consider the near horizon and near boundary 
form of the solution. First consider the near horizon behavior of solution 

<j){y) ^ -^L= Iniy - 1) + ^ ln(y - 1) (5.10) 
5V1 + C 2 

regularity at horizon implies that 



A = -^Vl + ^ (5-11) 

Take the limit y — y in fl5.8p and using (15. lip , we find the asymptotic form of 
solution near the boundary as 



2 2 

;y , y . y. 



</)(y) ~ 1 - iw{-b\l l + C'-+^-+^-) (5.12) 



and in original coordinate we have 



5i w 



EM \y^o- 1 + -r\P- + C'y^T y /yr (S-is) 
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